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A SIMPLE METHOD TO INCORPORATE WATER VAPOR ABSORPTION 
IN THE 15 pan REMOTE TEMPERATURE SOUNDING 

G. Dalu* 

C. Prabhakara 
B. J. Conrath 

ABSTRACT 

The water vapor absorption in the 15 pan C0 2 band, which 
can affect the remotely sensed temperatures near the surface, 
can be estimated with the help of an empirical method. This 
method is based on the differential absorption properties of 
the water vapor in the 11-13 pi m window region and does not 
require a detailed knowledge of the water vapor profile. With 
this approach Nimbus 4 IRIS radiance measurements are in- 
verted to obtain temperature profiles. These calculated pro- 
files agree with radiosonde data within about 2°C. 
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A SIMPLE METHOD TO INCORPORATE WATER VAPOR ABSORPTION 
IN THE 15 Mm REMOTE TEMPERATURE SOUNDING 

INTRODUCTION 

Temperature sounding, from remotely sensed infrared radiances in the 
:• i Mm C0 2 band, has been a significant advancement in the field of meteorology. 
Several investigators (Wark and Hilleary, 1969; Hanel and Conrath, 1969; Conrath, 
et al. , 1970; Smith, et al. , 1970; Wark, 1970) have developed techniques to in- 
vert the 15 Mm radiances to derive temperature profiles. In the 15 Mm band 
there is some absorption due to water vapor. In the earlier studies mentioned 
above, a water vapor profile was used to account for such absorption. 

Retrieval of the water vapor profile from remote measurements in one of 
the water vapor absorption bands, in turn, requires the temperature profile. 

Thus the solution to the problem of temperature sounding from the 15 Mm data 
becomes one of an iterative nature. It will be helpful to avoid this iteration. 

Water vapor decreases rapidly above the surface with a scale height of 
about 2km. Contrary to this fine scale, the vertical resolution of the retrieved 
temperature profile (Conrath, 1972) from the 15 Mm band is a few kilometers near 
the surface. In view of these facts, accurate representation of the water vapor 
profile near the surface is not necessary for the temperature sounding. Further- 
more the 11-13 Mm window measurements, which can be used ^ estimate the 
surface temperature (Prabhakars, et al. , 1974) allow us to estimate the 
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absorption due to water vapor. With this information, a simple method is de- 
veloped in the present study which utilizes the window measurements in the 
11-13 jum region to assess the water vapor absorption in the 14 (im wing of the 
CO 2 band without an interative procedure. The problem of temperature sounding 
is thus uncoupled from the water vapor sounding. 

Nimbus 4 IRIS spectra (Hanel, et al. , 1972) have been used to test the method 
developed, but comparable results are expected from multi-channel radiometer 
measurements. 

METHOD 

In Figure 1, a spectrum measured by IRIS over Guam on April 27, 1970 in 
the spectral region 600 to 1000 cm -1 is shown. 

Utilizing radiosonde measurements over Guam on that day, the absorption 
due to water vapor is calculated with a line by line integration program of Kunde 
and Maguire (1974) and is shown in the same figure. From this figure one can 
readily see the magnitude of water vapor absorption in the 15 Mm band which can 
be significant particularly over the tropics. Near the middle of the band, the 
C0 2 weighting functions peak in the stratosphere or upper troposphere (Conrat'j, 
1972). This can be easily inferred from the set of blackbody emission curves 
included in the figure. Hence the absorption due to the water vapor present in 
the lower troposphere is unimportant. However, the weighting functions for CO 2 
in the 14 Mm region of the band, from about 725 to 770 cm -1 , peak in the lower 
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Figure 1. Cloud-free Nimbus 4 IRIS spectrum taken over the ocean on 
April 27, 1970 at 15.1°N and 144. 7°W (GUAM), and the blackbody emis- 
sion curves. The water vapor transmittance calculated with the radio- 
sonde data is also shown. 
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troposphere or near the surface and the water vapor absorption should be taken 
into account. 

From the transmission curve of the water vapor displayed in Figure 1, one 
can also see that there is a simple trend to the curve from 11 to 14 pm. Since 
C0 2 has negligible absorption in the 11-13 pm window, radiance measurements 
in this region can give some information on the water vapor absorption in the 
atmosphere. Also the surface temperature can be estimated with these window 
measurements (Prabhakara, et al. , 1974, subsequently referred to as PDK). 
With the help of this information, we can assess the water vapor absorption In 
the 14 Pm region which is spectrally adjacent to the window. 

By considering the water vapor continuum absorption to be dependent on 
vapor pressure and total pressure, and also the absorption due to spectral lines, 
PDK have shown empirically that the transmission r in the window from 775 to 
960 cm’ 1 can be linearly related to the amount of precipitable water w(gcm -2 ) 
in the atmosphere as 

r„ = 1 -k„w (1) 

where k v (cm 2 g _1 ) is a relative absorption coefficient. Also from similar cal- 
culation of the water vapor transmission in the 725-775 cm -1 region, it is found 
that the approximation expressed by Equation (1) is valid. The values of k(v) 
used in this study are listed in Table 1. 
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In a cloud-free, non-scattering atmosphere under local thermodynamic 
equilibrium, the radiative equation for the upwelllng Intensity IO) may be 
written as 


I(»0 = BO', T 0 ) riv, p Q ) + 


f BO, 

•'riv, p ) 


T(p)l Mv, P) 


( 2 ) 


where p 0 and p are surface pressure and pressure at any height, respectively 
(mb), 

v is the wave number (cm -1 ), 

T is the temperature (K), 

T 0 is the surface temperature (K), 

B Ib the Planck intensity (erg ctn -1 ster" 1 s _1 ), 
r is the transmission from any pressure level p to the top of the 
atmosphere. 

The surface emissivity is assumed to be unity. 


For the window region, the integral in Equation (2) can be simplified and 
then the equation can be rearranged suitably to give 


where 


Bit', 


0) _ . . (. BQO V Bp) 

,T 0 ) ■ f0 ' ,Po] ) V ‘BO'.ToV B(*'.T 0 ) 


f Bp, 

*r{v, p 0 ) 

f p) 

Jr(v,p a ) 


T(p)l dri*'. p) 


(3) 


(4) 


k 
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j 

is the effective emission of the water vapor layer near the surface. 

In the window, the ratio of l(v )/B(v ,T 0 ) is a function of the total transmis- 
sion of water vapor in the atmosphere and the temperature distribution. The 
dependence on the temperature distribution is contained in B(i>). As the water 
vapor content in the lower atmosphere varies with height much more rapidly than 
the temperature, the water vapor distribution critically determines B(v). 

From Equations (1) and (3) we ha\e 

= (5) 

where 

Rw.-.V = »(l - <« 

The water vapor content w and the effective emission B(v) are not related in a 
simple fashion. However if we use the integration program of Kunde and Maguire 
(1974) for some typical average atmospheric models, we find that the water vapor 
content w and the function F(w, v, T Q ) are correlated in a simple way. In the 
graph of Figure 2 the values of w and F are plotted for the 831-887 cm -1 window 
region for an average tropical atmosphere (T on the graph), for an average mid- 
latitude atmosphere (MS) in the summer, and an average midlatitude atmosphere 
in the winter (MW). The temperature and water vapor profiles used for these 
model calculations are given in McClatchey, et al. , (1972). 
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Figure 2. Relationship between the function F and the precipitable water con- 
tent in the atmosphere. T is an example of tropical atmosphere, and MS and 
MW are midlatirude cases for summer and winter, respectively. 
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The surface temperature T 0 can be determined from a method developed by 
PDK with the help of two remote radiance measurements in the 11-13 Mm window. 
From the measured brightness ten .nature in the 831-887 cm -1 region, the ratio 
I(.v)/B(r, T Q ) can be calculated and with the help of Equation (5), the value of the 
function F(w, v , T ) is determined. The graph of Figure 2 then provides an eval- 
uation of the parameter w, which can be used as an equivalent water vapor content. 

Given the equivalent total water content w e , it is sumed that the water 
vapor w(p) contained from any pressure level p to the top of the atmosphere can 
be related to w e with the help of a scaling function R(p) such that 

w(p) = R(p! w c (7) 

The values used for H(p) deduced from some mean water vapor distribution are 
listed in Table 2. 

Table 2 

Scaling Function R(p) for the Water Vapor Distribution 


p(mbar) 

100 

200 

300 

400 

500 

600 

709 

800 

900 

1000 

I’(P) 

0.00 

0.01 

0.02 

0.05 

0.09 

0.15 

0.26 

0.44 

0.66 

1.0 


For the 15 Mm C0 2 band. Equation (2) can be written as 
\{u) = B(v.T 0 )t H 20 (p.p 0 >t ( t>2 <r.p 0 ) 


I. 


1 


B(r, T(p)l r co ^(v t P> dT n 2 o (,; ’ 


( 8 ) 


+ f B[*\ T( p > 1 r„ () (m, p»dr co (r, p) 

A-O^Po) 2 2 
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In this equation the intensity I(i') is measured; the value of B(»,T 0 ) is estimated 

by using at least two channels in the 11-13 Mm window region; the water vapor 

transmissivity r H ^ 0 (v, p) is estimated with the help of Equations (1) and (7); 

r ro O', p) is known, so the only unknown is T(p). Any known technique to in- 
2 

vert Equation (8) can be used to calculate the temperature profile. In the present 
study we have used the technique developed by Conrath (1972). 

In Equation (5), B is shown to be dependent on the temperature structure. 
However the equivalent water content estimated from Figure 2 can account for 
cases n which temperature gradients do not differ substantially from an average 
lapse rate of -6. 5°C/km. Only in the extreme cases where there is a superadi- 
abatic gradient or a strong inversion near the surface will this procedure be in 
error. Such a limitation is also inherent in the more elaborate iterative inver- 
sion techniques, because of the crude vertical resolution of the retrieved tem- 
perature profiles (Conrath, 1972). 

RESULTS 

The Nimbus 4 IRIS spectra were used to test the method developed. A few 
spectra have been selected for which the corresponding temperature profiles 
from radiosondes were available. The radiosonde launches were conducted 
within one hour and a few degrees in latitude and longitude from the spot sensed 
by the satellite. 
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As an example, the temperature profiles obtained near GUAM, latitude 15°N 
and longitude 145°E, on April 27, 1970, are presented in Figure 3. The solid 
line corresponds to the radiosonde temperature profile, while the dashed line is 
the corrected temperature profile obtained from the inversion of the 15 pm C0 2 
band. If no correction for the water vapor absorption is introduced, i.e. , this 
absorption is assumed to be negligible at any level, the inversion of the C0 2 band 
would give the deformed temperature profile shown with the dots in the same 
figure. 

In Table 3 a few temperature profiles are listed. For each case there are 
three profiles: the first one was obtained with a radiosonde; the second one was 
calculated from the 15 Pm C0 2 band, but corrected using the water vapor profile 
obtained from the radiosonde itself; the third one was calculated with the present 
method, i.e. , without any data on the water vapor content of the atmosphere. 

The error of the temperatures calculated from the inversion of the radiances 
measured from the IPIS, is about 2°K in the layers near the surface. 

The correction for the water vapor absorption is actually most important 
between TOOmbar and the surface (Figure 3). The comparison of the profiles 
presented in Table 3 show differences of 1°K or less between this method and 
the inversion of the 15 Pm band with the correction f ? the water vapor absorp- 
tion obtained by using the water vapor profile measured with a radiosonde. 
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Figure 3. Comparison among the radiosonde temperature profile, the retrieved 
profile without water vapor correction, and the corrected temperature profile. 
The IRIS data used are the same presented in Figure 1 (GUAM). 
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Comparison of the Present Temperature Sensing Method with the Radiosonde Measurements 
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CONCLUSION 

The method developed to acoount for the water vapor absorption in the in- 
version of the 15 Mm C0 2 band has given satisfactory results. It has been found 
that the temperature profiles calculated with this method are very close to the 
temperature profiles obtained using the measured water vapor profile. 

The main advantage of the method proposed is that the same two channels in 
the 11-13 Mm window region needed to obtain the corrected surface temperature 
can be used in a simple scheme to oorrect the temperature profile for the water 
vapor absorption effect. 

The calculated temperature profiles agree with the radiosonde temperature 
profiles within the accuracy of the measurements, and show an even better agree- 
ment with the profile corrected with the measured water vapor profile. 

According with the results obtained, a nine channel radiometer, incorporat- 
ing seven in the 15pm C0 2 band and two in the 11-13 Mm water vapor window 
region, can provide a satisfactorily accurate remote temperature sounding. 
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FIGURE CAPTIONS 

Figure 1, Cloud-free Nimbus 4 IRIS spectrum taken over the ocean on April 27, 
1970 at 15. 1°N and 144. 7°W (GUAM), and the blackbody emission curves. 
The water vapor transmittance calculated with the radiosonde data is also 
shown. 

Figure 2. Relationship between the function F and the predpltable water con- 
tent in the atmosphere. T is an ej ample of tropical atmosphere, and MS 
and MW are midlatitude cases for summer and winter, respectively. 

Figure 3. Comparison among the radiosonde temperature profile, the retrieved 
profile without water vapor correction, and the corrected temperature pro- 
file. The IRIS data used are the same presented in Figure 1 (GUAM). 



